NOTICE 


THIS DOCUMENT HAS BEEN REPRODUCED FROM 
MICROFICHE. ALTHOUGH IT IS RECOGNIZED THAT 
CERTAIN PORTIONS ARE ILLEGIBLE, IT IS BEING RELEASED 
IN THE INTEREST OF MAKING AVAILABLE AS MUCH 
INFORMATION AS POSSIBLE 



DOE/NASA/1 0769-1 5 
NASA TM-81 756 


Preliminary Investigation of 
Acoustic Oscillations in an H2-O2 
Fired Hall Generator 


(MJkSA'TH-81756) PBBLIHIHABI 18VBST1GAT10H H81-23610 

OF ACOUSTIC OSCILLATIOMS IB AB H2>02 FIBBD 
HALL GBBBBATOB (NASA) 13 p HC A02/HF A01 

CSCL 10A Onclas 

G3/44 42352 


Bert Phillips 

National Aeronautics and Space Administration 
Lewis Research Center 


Work performed for 

U.S. DEPARTMENT OF ENERGY 
Fossil Energy 

Office of Magnetohydrodynamics 


Prepared for 

Nineteenth Symposium on the Engineering 
Aspects of Magnetohydrodynamics 
Tullahoma, Tennessee, June 15-17, 1981 




DOE/NASA/10769-15 
NASA TM-81 756 


Preliminary investigation of 
Acoustic Osciiiations in an H2-O2 
Fired Haii Generator 


Bert Phillips 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 


Performed for 

U S. DEPARTMENT OF ENERGY 
Fossil Energy 

Office of Magnetohydrodynamics 
Washington, D.C. 20545 

Under Interagency Agreement DE-AI01-77ET10769 


Nineteenth Symposium on the Engineering 
Aspects of Magnetohydrodynamicj. 
Tuliahoma, Tennessee, June 15-17, 1981 


PftELIMINMY INVESTIGATION OF ACIKISTIC OSCILLATIONS 
IN AN H2-02 FIRED NALL GENERATOR 

Nrt Rhilllps 

NASA Lewis Rtstarch Center 
Clevelend, Ohio 4413S 


Abstract 

In this paper are presented burner pressure 
oscillations and Interelectrode voltage oscilla- 
tions Measured In an open-cycle supersonic flow 
Hall generator. The lonlaed gas for the channel 
was supplied by seeding the approxinately 1 Ib/sec 
of hydrogen-oxygen coMbustlon products with cesi- 
um. Since both the burner and the channel were 
located within Magnetic fields exceeding Tesla 
during operation, an Infinite probe pressure Mea- 
surement technique was used to measure burner 
pressu.'e oscillations. Calibration of the burner 
pressurr transducer using a resonance tube tech- 
nique Is presented. Evidence is presented for the 
existence of the first longitudinal Mode of os- 
cillations (SOOO Hi) within the burner. Inter- 
electrode voltage oscillations were simultaneously 
measured at two separate axial stations. The 
magnitude change and the phase shift between the 
two signals was Interpreted as a decaying magneto- 
acoustic wave driven by the burner that propagates 
at local gas plus sonic velocities. The amplitude 
of the electrical voltage oscillations at the 
start of the power producing region of the channel 
varied with the magnetic field. This variation Is 
compared with the results of a simple perturbation 
analysis. Arguments are presented for using an 
unsteady model for analyzing wave processes In 
channels. 

1. Introduction 

An Important property of an electrical genera- 
tor Is Its stability. In a formal sense, the 
stability can be recognized as the sensitivity of 
the generator to disturbances that will cause a 
time varying output In the electrical power. An 
MhD generator with Its closely coupled fluid mech- 
anical and electrical fields Is particularly sus- 
ceptible to a variety of Instabilities. These can 
arise from variations In operation of one or more 
of the systems connected to the generator, or they 
can arise from the generator Itself. In order to 
achieve a power generation capability suitable for 
commercialization. It ’s Important that these In- 
stabilities be understood and their effects mini- 
mized. Open-cycle ftiD generators are particularly 
sensitive to disturbances that arise within the 
burner that supplies the Ionized gas to the gen- 
erator. Depending on whether the channel is 
operated as a sub- or supersonic generator, any 
burner-generated disturbances may be able to cou- 
ple back to the burner, resulting In sericus sys- 
tem Instabilities with potentially damaging re- 
sults. While there has been no evidence presented 
in the literature for the existence of these dam- 
aging Instabilities, It Is useful tv recall the 
history of development of the la>gc liquid and 
solid propellant rocket engines used for space 
propulsion. The smaller engines, that were used 
for preliminary and subscale development, were 
either stable or the Instabilities were not seri- 
ous. Only when the larger engines were under 


development did we see the destruction that high 
aMplltude pressure oscillations could cause. The 
Increased iMportance of these oscillations as the 
scale Is Increased Is related to the decrease in 
the relative wall damping and the proximity of the 
burner or engine time scales with those of the 
combustion processes. 

In order to experimentally study the nature of 
Instabilities In open-cycle NhD generators. It Is 
convenient Initially to decompose what can be a 
very complex problem Into Its more simple compo- 
nents. For the case of burner-coupled oscilla- 
tions, this car be achieved most readily by oper- 
ating the generator supersonically. Thus any dis- 
turbance In the burner can be used to study the 
ability of the generator to aiqillfy or decay It 
without a feedback to the burner. The most obvi- 
ous requirement for this study Is a source of dis- 
turbances within the burner. Either a pulse or 
oscillatory disturbance would be suitable. Based 
on the requirements Stated above. It appeared that 
the H 2 -O 2 Fired (*iD Generator Experiment at 
the NASA Lewis Research Center would be suitable 
for this study. The output of the generator was 
of the order of 100 MW/mJ, the channel was cap- 
able of Instrumentation for transient measure- 
ments, and, most Important of all, there was son« 
evidence that the natural frequency for the first 
longitudinal mode of oscillation with the combus- 
tor was in tiie range of a few kilocycles which 
would be easily measured. Therefore, this mode 
could be used as a source of oscillation propagat- 
ing into the channel. 

11. Hq -Q.. Burner and Channel 

The mo channel used In tnis Investigation has 
been more completely described In previous 
papers. 1 Briefly, It consisted of a gaseous 
hydrogen, gaseous oxygen fueled water-cooled 
rocket combustor, operating stolchlometrical ly, 
with a nominal flow rate of 1 Ib/sec, a chamber 
pressure of 10 atmospheres, and seeded with a 
cesium hydroxide water solution. Both steady and 
transient combustor pressures were measured. The 
nozzle exit (channel entrance) Mach number is 2 
with a diameter of 2 Inches. The channel is a 
segmented Mall design with half-inch uncooled cop- 
per heat sink electrode plates separated by mica 
Insulators. The channel Is approximately 2 feet 
long with a conical bore, a 4 to 1 area ratio, and 
a single resistance load. Steady-state static 
pressures and voltages are measured at various 
axial positions down the channel. Oscillatory 
voltage signals were measured at two positions. 

The first position was at the start of the pox- 
producing region of the channel and the other at a 
position ~6 Inches downstream. The liquid neon 
cooled cryomagnet used can achieve B fields In 
excess of 5 Tesla for short (5 sec) runs. The 
maximum field strength used In tns Investigation 
was 4 Tesla. A pictorial sxetch of the magnet, 
generator, and burner used is shown In Fig. 1. 


111. Weisurement of Cnamber 
Pi^essure Oscniatiofis 

In order to characterize the pressure oscllli- 
tions within a comoustor, it is necessary. In gen- 
eral, to have one or more high frequency pressure 
transducers flush mounted in the changer wall. 
Thus, for normal rocket engine operation, the 
transducers must be cooled, and if necessary, 
mounted to avoid the effects of wall vibrations on 
the signal. In the case of the channel being 
s!ii:iied, the channel together with the burner is 
irserted within the narrow magnet bore. Both the 
high magnetic fields present during operation and 
space limitations prevented installation of such a 
flush mounted transducer. Therefore, in order to 
measure the chamber pressure oscillations. It was 
necessary to separate the transducer from the 
chamber by a long small diameter tube. A review 
of the literature indicated that the most promis- 
ing approach was to use a very long tube with the 
high frequency pressure transducer flush mounted 
to the inside of the tube ("infinite probe") as 
shown in Fig. 2. In principle, the damping 
created by the tubes excess length would prevent 
establishing any starring wave patterns and give a 
flat frequency response curve up to the poirt 
where viscous and thermal dissipation would cause 
significant signal decay. In Ref. 2, the authors 
present a technique for including damping material 
within the tube at a location calculated to elimi- 
nate any standing - patterns and insure a flat 
response. Althoug .his technique was demon- 
strated to be suitable for the oscil' ’ons stud- 
ied in that reference, it was not pos lole to 
include damping material in the present study 
since the tube was used as part of the steady- 
state pressure measurement and burner purge sys- 
tem. Consequently, it was necessary to provide a 
means of calibratirc the pressure transducer in 
its "infinite probe" configuration. 

Burner Pressure Calibration 

A convenient method for the generation of high 
frequency pressure oscillations together with ele- 
vated mean pressures is the resonance tube (Fig. 

This Simple device depends on the initiation 
of longitudinal organ pipe oscillations within a 
'■nail dosed-ended cylindrical cavity by having a 
high pressure gas jet cinted into the open end of 
the cavity at specifi; distances fruii the open 
e'lj. This device has been extensively studied and 
1 ' well characterized. For a frequency range of 
Kuo to 5000 HZ, the anticipated frequency range 
for any burner pressure oscillation, cavity depths 
of \I2 to 3 inches are required, using air as the 
driving gas. 

Frr a selected pressure of the air, the spac- 
ing batween the air jet and the cavity inlet is 
pres'/t at a point where maximum oscillations 
occi.r. Frequency of oscillations is varied by 
adjusting the position of the cavity piston. 

With a resonance tube diameter of 1/4", it was 
possible to achieve oscillation frequencies in 
excess of 6500 Hz using 125 psig air as the driv- 
ing gas. In order to assure that the oscillations 
generated oid not have a significant component of 
harmonic distortion, which would have complicated 
the calibration process, the output from the flush 
mounted pressure transducer was displayed on an 
oscilloscope screen and photographs were taken of 
the wave shape. A typical example is shown in 
Fig. 3. The Shape is remarkably sinusoidal con- 


sidering that the anplltudc of the signal 1y 15 
psi peak-to-peak out of a mean pressure of 40 
psia. The frequency of the displayed wave Is SOCiO 
HZ with some small variation about the value as 
seen by the thickness of the trace. The amplitude 
of the pressure oscillations varied from 10 to 15 
psi as the frequency of the wave was Increased 
from 1000 to over 6500 Hz, providing a relatively 
constant input for calibration purposes. In addi- 
tion, the static pressure within the tube, as mea- 
sured by a Bourdon gauge at the end of the In- 
finite probe coll, was 40 psia. 

The signals from the pressure transducer flush 
mounted to the inside of the resonance tube and 
the transducer attached to the infinite probe were 
recorded on magnetic tape. 

The recorded signals were then analyzed using 
a Hodel 400B Nicolet Scientific Corporation Fast 
Fourier Transform Analyzer. This digital-type 
analyzer readily produces the time variation of 
the individual signals Power Spectral Density as 
well as the Cross Spectral Density of the two sig- 
nals. The result Is a capability to observe the 
time variation of the phase angle between the two 
signals, the ratio of the magnitude of the two 
signals or transfer function, and a measure of 
the correlation between the two signals, i.e., the 
coherence. The transfer function and the coher- 
ence for the recorded data are shown in Figs. 4 
and 5. Up to 4000 Hz, the transfer function indi- 
cates a series of resonant peaks within the probe, 
pernaps due to wave reflections off of area dis- 
continuities, superimposed on a straight line 
decrease due to viscous and thermal losses within 
the probe. Beyond 4000 Hz there is a break in the 
slope of the viscous/ thermal loss curve. Referr- 
ing now to the coherence l>etween the two signals, 
shown on Fig. 5, it can be seen that the coherence 
drops off very rapidly beyond 4000 Hz, indicating 
that data in this area are either insufficient or 
that the signals are uncorrelated. Rather than 
repeat the calibration process using the digital 
analyzer, both signals were displayed on an os- 
cilloscope screen and a series of photographs were 
taken of the wave shapes. The amplitude ratio of 
the waves, based on the oscilloscope photographs, 
is given as t c circular symbols on Fig. 4. At 
frequencies less than 4000 hz, the two methods of 
analysis agree whereas at frequencies above 4000 
HZ, the scope photographs confirm a constant slope 
in the transfer function due to viscous and ther- 
mal losses in the infinite probe. Based on Fig. 

4, an Infinite probe transfer function of 28 d§ 
was used to convert the measured 5000 Hz chamber 
pressure oscillations to actual values. 

Burner Pressure Oscillations 

During a test, the output from the infinite 
probe transducer was recorded on magnetic tape for 
analysis of the pressure oscillations. Using the 
analyzer previously discussed, plots of power 
spectral density as a function of either time or 
averaged over the 5 second run could be readily 
obtained. A typical time averaged plot is shown 
in Fig. 6. The magnitude of the peak at approxi- 
mately 5000 HZ is several dB above the background 
and corresponds to a measured pressure oscillation 
of 0.84 psi (p.t.p.), providing an unmistakable 
signature of burner combustion oscillation. The 
burner was cylindrical in cross section, 2.5 
inches in diameter with a distance from the injec- 
tor to the throat of seven inches. Based on a 
calculated chamber Mach number of 0.23 and adopt- 
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Ing tht convention that the pressure wave will 
reflect from a position n1dwa> between the start 
of the convergent section of the burner chamber 
and the throat, a first lonoltudlnal mode resonant 
frequency of 5262 Hz was calculated. Considering 
the assumption of the proper reflection point for 
an axial pressure wave In a converMnt duct, the 
agreement between the experimentally measured fre- 
quency and that predicted confirms the nature of 
the oscillation mode within the burner. 

The amplitude of the oscillations, using the 
previously discussed calibration procedure, was 
21.3 psi peak-to-peak or approximately 15 percent 
of the mean chmnber pressure of 10 atmospheres. 

IV. Measurement of Channel 
Electrical dsci11at1ons ~ 

Hall Voltage Oscillati ons 

During a typical test run, the vu’tage between 
two sets of electrodes one set at the start 
the power producing region of the channel and th,' 
other approximately 6 Inches downstream, was re- 
corded on magnetic tape. Using the digital ana- 
lyzer previously discussed, a variety of ap- 
proaches were available to analyze the data. 

Power spectral density plots were used to Identify 
the oscillatloi. frequency and determine the 
signal-to-noise ratio. A typical plot is shown In 
Pig. 7 for the two recorded signals. While there 
Is considerable noise, the peaks at 5000 hz are 
readily distinguishable. In this case, the 
signal-to-noise level is of the order of 3 dB 
only. The effect of the high noise level is to 
decrease the coherence of the signals for phase 
lag measurements. Both coherence and phase lag 
are shown in Fig. 8. For the frequency of inter- 
est, the coherence is in excess of 0.50. The 
averaged value of the phase lag for that frequency 
was taken as 72*. 

Estimate of Channel Wave Propagation Velocity 

In order to determine whether the oscillations 
from the combustion chamber propagated at the 
local gas velocity in the channel, corresponding 
to entropy wave production within the burner, or 
at the sum of both gas and local sonic velocity, 
corresponding to an acoustic wave. It is necessary 
to compare the measured wave velocity based on the 
electrical signal phase lag, the frequency, and 
the distance between the measurement stations with 
an estimate of the gas and sonic velocities within 
the channel. The measured wave velocity based on 
phase angle and frequency was, from Table I, 

11,250 ft/ sec. 

A one-dimensional analysis of channel perfor- 
mance was used to estimate the channel conditions 
that could not be readily measured. The results 
of the channel analysis are shown In Fig. 9 for a 
nominal B field of 3 Tesla. On the plot are shown 
the variations of static pressure. Hall voltage, 
and the sum of gas plus sonic velocities. For the 
acoustic wave velocity, pressure, and steady-state 
Hall voltage, sufficient data were taken to assure 
the validity of the calculation procedure. The 
average wave velocity between the two measurement 
points as indicated by the figure and Table I Is 
11,075 ft/sec which agrees sufficiently well with 
the theoretical number to verify that the elec- 
trical wave propagates as an acoustic-like wave 
rather than an entropy wave. 


Relative Magnitude of the Electrical OKlllatlows 

FKlllty limitations prevented the measurement 
of oscillatory pressures et the same position as 
the electrical signal so that In order to deter- 
mine whether the data were meaningful. It was 
necessary to try to simulate the electrical signal 
amplltute. The first or upstream position at the 
front of the power producing region wts 6.6 Inches 
from the channel entrance and correspond to 
0.5 Inches downstream from the electrodes used for 
current lead-ins. Based on the short distance and 
without any other Information, a simplified per- 
turbation model with no wave damping mechanisms, 
based on the work of Narble and Cande1,A was 
used to estimate ttM osc111a*1ons In pressure, 
velocity, and temperature a. that location, given 
only the magnitude of the chamber pressure os- 
cillations. Using the value from the burner pres- 
sure measurement, a relative magnitude (oscillat- 
ing Hall voltage/ steady Hall voltage) of 0.258 was 
calculated. This value, as shown on Table I com- 
pares very favorably with the experimentally mea- 
sured value of 0.268. The agreement Is probably 
fortuitous, considering lEhe assumptions Involved, 
but does Indicate that the electrical oscillations 
near the entrance to the channel are explainable. 

Variation of the oscillations at the first 
axial position with ft.poscd B-fleld was also cal- 
culated using the perturbation analysis already 
discussed. A con^arlson of the analyses with 
experimental data Is shown In Fig. 10. Both the 
data and simple perturbation theory Indicate that 
the relative level of the oscillations at that 
position decrease as B field Is Increased. The 
one-dimensional channel model incorporated in the 
perturbation analysis Is not valid for very low B 
fields so that the curve was not continued. Addi- 
tional measurements and analyses are necessary to 
verify this trend. 

Oscillatory Pressure and Electrical Signal Decay 

The results shown In Table I for the change In 
relative magnitude (oscillating Hall voU- 
age/ steady Hall voltage) of the electrical os- 
cillation with axial position in the channel are 
related to the loss mechanisms for oscillations 
within the WiD duct. Oscillatory loss mechanisms 
are generally based on heat transfer and friction- 
al losses CO the walls together with Kh0-r«1ated 
amplification or loss mechanisms. One of the 
approaches commonly taken^ to analyze these 
phenomena Is to assume that the frictional and 
heat loss mechanisms are the same for the oscilla- 
tory flows as for the steady flows (quasi-steady 
assumption). The validity of the assumption can 
be viewed In terms of the ability of the momentum 
thermal boundary layers to respond to osHlatory 
changes. Boundary layers have a finite time to 
react to variations In the external potential flow 
that Is proportional to the square of the boundary 
layer (B.L.) tnickness divided by the kinematic 
viscosity. The proportionality factor Is of the 
order 0.1, based on exact solutions of boundary 
layer startup and boundary layer oscillations. 

The ratio of B.L. relaxation time to wave time 
should then provide a measure of the validity of 
the quasi-steady assumption for analyzing wave 
dampino processes. For the case of the present 
channel study, the ratio of B.L. relaxation 
time/wave time Is approximately 50:1. Thus the 
quasi-steady assumption does not appear to be 
valid for analyzing this system. As a comparison, 
a 500 Mb thermal coal-fired channel with a 50 msec 
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rtsttftnce t<ne ceal'firctf cMbuytor Mill bm • 

I.L. time wm t)m« ratio of tbo orOtr of 100 to 
300 doptnOiot on tbo natoro of tbo okI nations 
and otbor factors. Thus, It apfoars to M ntctt> 
ury to dovolop osttMtos of tbo unstoa^y lA. 
lots mtcbanUirs In M’Otr to accoratoly cbarac- 
tcrlio wave procotsot In tboso cbwMois. 

¥, OlscitattOB 

Tbo rosuitt ^talnod from tbis study Indicata 
that tbo acoustic waves prodiicod by a caobuttor 
^r«tinp at the first loAOltndtiia) ando of ot- 
cl nation ^ay within the duct. Willo the aeipll- 
tubes of tnt oscillations In Halt vottapa at the 
first or upstream statim art consistont with a 
simplified perturbation a«dtl. analysis of the 
decay of the otc it tat ions down the channel should 
be based on an unsteady bonndwy layor amdtl 
rather than a quasi-steady madol. In addition, 
the use of an Infinite probe aponoach for measur- 
ing pressure oscillations within an MtO c omb u stor 
has been assessed as being reasonable provided 
that seme means for calibrating t‘« pi^ can be 
found. The resonance tube tfipnteh for calibra- 
tion ^pears to be suirsble for this purpose. 
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TablO 1. - NAU VOLTABE AMALVSIS 
Cl fltld ■ 3 Tesla, f re q uenc y of oKltle- 


tion ■ iobo Mt amplitude v burner pres- 
sure Mclilatloh • 21 psl Cptp).] 

location of transient 
eeltoge sensor 
(Inches frem inlet) 


12.0 

^eedy-state 

Imerelecteede 

uoltege 

20. 

32. 

Transient 

Intereleetrede 

voltage 

S.37 

3.25 

Tr«it1entfsteady 
ratio (aaperlmental) 

.268 

.101 

Trans lent/steedy 
ratio (theoretical) 

.258 


Phase lag 

7Z degrees (0.00004 sec) 

Nave velocity 
(experimental) 

11,250 ft/sec 


Have velocity 
(theoretical) 

11,075 ft/sec 
(Mch No. > 1.86) 
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Figure 1 . - GH2-GO2 combustion MHD experiment installation. 



Figure Z - Resonance tube calibration apparatus for infinite probe. 
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Figure 4 - Infinite prabe transfer function. 
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Figure Id - Ratio of h^nsient/steady interelectrode voltage 
(upstream posKkm). 


